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Abstract.  The  influence  of  collision  broadening  due  to  phonons  and  the  intra- 
collisional  field  effect  on  the  inverse  Auger  process  have  been  examined.  Results 
are  presented  for  silicon  and  show  for  the  first  time  the  wave-vector  dependence 
of  the  ionization  rate  including  collision  broadening. 


1.  Introduction 

Impact  ionization  (the  exact  inverse  of  the  Auger  process) 
in  semiconductors  has  been  studied  extensively  since  the 
first  theoretical  investigations  of  hot  carrier  transport.  An 
accurate  description  of  the  ionization  process  in  solids 
requires  a  full  band  structure  treatment  along  with 
associated  constitutive  complexities;  inclusion  of  band 
structure  effects  requires  detailed  numerical  calculations 
to  establish  the  magnitude  of  the  matrix  elements  for  the 
process.  Initial  efforts  to  understand  impact  ionization 
did  not  consider  detailed  band  structure  effects,  but 
developed  simple  phenomenological  models.  Prominent 
among  initial  studies  of  impact  ionization  was  the 
Keldysh  model  [1],  originally  developed  for  direct  gap, 
parabolic  band  materials;  this  model  has  been  used 
extensively  in  applications. 

Unfortunately,  the  Keldysh  formula  is  not  applicable 
for  most  real  materials  due  to  the  complexities  intro¬ 
duced  by  band  structure  effects;  specifically,  the  formula 
does  not  adequately  describe  indirect  gap  materials  such 
as  silicon.  Kane  [2]  has  performed  a  more  exact,  Fermi 
Golden  Rule  calculation  of  the  ionization  process  includ¬ 
ing  a  realistic  band  structure.  However,  the  Fermi 
Golden  Rule  does  not  adequately  describe  the  impact 
ionization  process  for  two  reasons;  first,  because  the 
phonon  scattering  rate  is  typically  very  large  for  the  high 
energy  electrons  which  undergo  impact  ionization,  the 
lifetime  of  the  electron  states  is  short,  thereby  broadening 
the  ionization  transition;  secondly,  electric  fields  are 
typically  very  large  so  that  the  ionizing  electron  can  be 
accelerated  considerably  during  the  period  of  collision, 
thereby  affecting  the  dynamics  of  the  collision  process. 

In  this  paper  a  theory  of  impact  ionization  in  semi¬ 
conductors  is  presented  which  extends  the  initial  work  of 
Kane  by  including  the  influence  of  high  electric  fields  and 
high  scattering  rates  on  the  electron -electron  collision 
process,  and  utilizes  a  more  advanced  treatment  of  band 
structure  and  dielectric  response  screening  than  was 
available  to  Kane  in  his  original  work  [3].  Hence,  the  use 
of  Monte  Carlo  methods  combined  with  density  matrix 


formalism  results  in  a  major  extension  of  Kane’s  theory 
beyond  the  Golden  Rule.  In  particular,  this  investigation 
shows  that  the  already  soft  threshold  of  the  ionization 
rate  as  indicated  in  Kane’s  theory  is  considerably  broad¬ 
ened  by  electron-phonon  collisions  and  the  intra- 
collisional  field  effect,  so  that  a  well  defined  threshold 
does  not  exist.  Furthermore,  results  show,  for  the  first 
time,  the  wave  vector  anisotropy  of  the  ionization  rate 
in  silicon;  the  importance  of  wave  vector  dependent 
thresholds  is  also  discussed. 


2.  Quasi-particle  Boltzmann  equation  for 
impact  ionizing  electrons 


Many  authors  have  attacked  the  problem  of  finding 
quantum  transport  equations  which  improve  upon  the 
Boltzmann  equation  in  the  high  field  limit  [4].  In  this 
section,  we  briefly  outline  a  theory  which  includes  the 
effects  of  the  electric  field  and  phonon  scattering  on 
impact  ionization.  A  more  detailed  description  is  given  in 
[3].  In  order  to  incorporate  the  electric  field,  we  choose  a 
gauge  based  on  the  vector  potential,  namely  (see  for 
instance  [5, 6]), 


-4(f)  =  -c 


F(t’)  df' 

0 


ID 


where  F  is  the  electric  field  and  <  is  the  speed  of  light.  We 
can  label  our  states  with  the  usual  Bloch  band  indices 
and  accelerated  Bloch  wave-vectors  K  such  that 


K(k.  0  =  k  +  Eyt  Fv  =  (e  h)F  (2) 

where  k  is  the  usual  Bloch  wave-vector.  This  trans¬ 
formation  allows  an  exact  inclusion  of  the  electric  field. 
The  Hamiltonian  is  written  as 

W,o,  =  Wo  +  W  b  +  War  +  V  (3) 

where  H0  contains  the  crystal  potential  and  the 
electronic  kinetic  energy,  is  the  phonon  bath  Hamil¬ 
tonian,  Hab  is  the  system  bath  interaction,  and  Fis  the 


man,  10.10  ;QO  /AORKOfi  -i.  A3  tna  KA  rd\  1092  IQP  Publishing  Ltd 
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electron-electron  interaction  corresponding  to  impact 
ionization. 

Then,  we  look  for  an  equation  of  motion  for  the 
diagonal  of  the  electron  density  matrix  which  corre¬ 
sponds  to  the  semi-classical  distribution  function.  This  is 
accomplished  by  applying  the  following  projection 
operators  (see  for  instance  [7]) 

p  =  ^d®sK-.*  Q  =  i-r  (4) 

irBc 

to  the  equation  of  motion  for  the  density  matrix, 

ih^  =  [H,oMp]  (5) 

and  back-transforming  the  result  from  the  accelerated 
state  basis  to  the  stationary  Bloch  basis. 

The  resulting  equation  gives  an  expression  for  the 
impact  ionization  scattering  rate  R(  1 2  — ►  34,  r)  from  an 
initial  state  \nlkl;n2k2y  to  a  final  state  |n3lc3;  n4fc4>, 
which  includes  virtual  phonon  processes  and  accelera¬ 
tion  in  the  electric  field.  Note:  for  large 
r,  R(12-*34,  t)  -» /?( 1 2  -*  34).  This  exression  replaces 
the  energy  conserving  delta  function  of  the  Fermi  Golden 
Rule: 

R(12-34,f)  *  J'dr'S(I2.  (6) 


S(12,  34;  i.  /’)  =  exp 
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Figure  1 .  Impact  ionization  rates  /?av„(£)  for  Si  averaged 
over  initial  electron  energy  measured  from  the  bottom  of 
the  conduction  band.  Solid  curve,  zero  field,  no  collision 
broadening;  dotted  curve,  collision  broadening  (CB); 
chain  curve,  F  =  5  x  1 0®  V  cm  ’  and  collision  broadening. 


Here  r(nf,  kt)  and  A(n,,  kt)  are  the  imaginary  and  real 
parts  of  the  electron  self-energy  (due  to  virtual  phonons), 
and  |  F]2;34|  is  the  matrix  element  for  the  electron -electron 
interaction.  The  total  scattering  rate  for  an  initial 
electron  R{nlk1)  is  given  by 

*(«i*i)  =  1*02-34)  (8) 

34 

and  the  rate  averaged  over  initial  electron  energy, 
*ave(£k  >s 

P  _  E.jBzd3*,  *(«,*,)<$(£„(&,)-  £) 

Mazd  %5(EM-E)  ■  W 

The  expressions  above  are  used  to  calculate  the 
impact  ionization  rate  for  silicon  under  different  field 
conditions. 


3.  Quantum  effects  and  ionization  thresholds 

First,  we  examine  the  threshold  dependence  of  Rave(£)  for 
silicon.  Figure  1  shows  Rave(£)  calculated  using  the 
Fermi  Golden  Rule,  and  deviations  from  the  Fermi 
Golden  Rule  which  occur  because  of  the  inclusion  of 
either  collision  broadening  alone,  or  from  collision 
broadening  and  the  intra-collisional  field  effect  operative 
simultaneously;  thus,  deviations  from  the  Fermi  Golden 
Rule  results  are  due  to  the  inclusion  of  specific  time 
dependent  processes  of  finite  time  duration.  The  Golden 
Rule  result  shows  a  rapid,  threshold-like  increase  in  the 
scattering  rate  just  below  1.5  eV;  hence,  £lh  %  1.5eV. 
However,  when  collision  broadening  is  included,  the 
‘threshold’  decreases,  allowing  impact  ionization  to  occur 
well  below  the  silicon  band  gap.  This  clearly  shows  the 
influence  of  energy  broadening  due  to  the  phonons.  The 
ionization  rate  remains  finite  for  energies  down  to  0.8  eV. 
Lower  values  of  initial  energy  could  not  be  computed 
because  of  numerical  difficulties.  For  higher  electron 
energies,  the  effect  of  collision  broadening  is  relatively 
small;  collision  broadening  has  its  greatest  influence  near 
thresholds  where  the  rate  increases  rapidly  with  energy. 

In  the  presence  of  high  fields  (£  =  105  Vcm~ l),  the 
lower  end  of  the  ionization  curve  is  broadened  even  more 
dramatically.  Depending  on  the  orientation  of  the 
electric  field  with  respect  to  band  structure,  an  impact 
ionizing  electron  can  either  gain  or  lose  energy  from  the 
field  during  the  period  of  a  collision.  The  result  is  that 
low  energy  electrons,  unable  to  impact  ionize  in  the  zero- 
field  limit,  can  gain  enough  energy  from  the  field  during 
the  time  of  the  collision  to  ionize;  this  process  can  be 
viewed  as  an  effective  two-particle  Franz- Keldysh  effect. 
The  effect  is  again  most  dramatic  for  electrons  near  the 
Golden  Rule  threshold. 


4.  Wave-vector  dependent  thresholds 

Anderson  and  Crowell  [8]  were  the  first  to  consider  the 
influence  of  crystallographic  direction  on  threshold 
energy;  an  electron  with  a  given  initial  wave-vector,  kt. 
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Figure  2.  (a)  Ionization  rate  in  the  <111)  direction. 
Upper  set  of  curves:  second  conduction  band.  Lower  set: 
first  conduction  band.  Solid  lines,  collision  broadening 
R(k);  dotted,  no  collision  broadening  R(k):  chain, 
ffave (£(*))•  (b)  Energy  in  <1 1 1  >  direction;  solid,  second 
conduction  band;  dashed,  first  conduction  band. 
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Figure  3.  (a)  Equal  ionization  rate  curves  in  the  kz  =  Q 
plane,  first  conduction  band,  (b)  Equal  energy  curves  in 
the  k,  =  0  plane,  first  conduction  band.  The  energy  curves 
are  in  increments  of  I.OeV.  The  ionization  rate  contours 
correspond  to  the  average  rate  for  the  equal  energy 
surfaces  flave(1  eV),  /?„,„( 2eV),  etc. 


can  impact  ionize  only  if  its  initial  energy  exceeds  a 
minimum  energy,  £lh(k,),  where  fc;  is  the  unit  vector  in  the 
k,  direction.  Computations  of  £lh(k,)  in  references  [8] 
and  [9]  indicated  a  highly  anisotropic  threshold  for  Si  of 
about  3  eV  in  the  <  1 1 1  >  direction  while  in  the  <  1 10)  and 
<100>  directions  thresholds  of  2.1  eV  and  1.1  eV  were 
found.  This  means  that  rates  with  explicit  k  dependence, 
R(k),  could  deviate  significantly  from  Rave(£(k)). 
However,  the  above  referenced  calculations  have  as¬ 
sumed  that  the  wave-vectors  of  the  particles  involved  in 
threshold  processes  are  all  parallel  to  the  initial  particle 
wave-vector.  This  assumption  is  too  restrictive,  and 
overestimates  the  wave-vector  anisotropy  significantly. 

In  this  paper,  we  have  calculated  R(k)  on  a  0.05  x  2n/a 
mesh  in  silicon,  with  and  without  collision  broadening; 
the  results  of  this  calculation  are  shown  in  figures  2  and  3. 
Figure  2(a)  shows  R(fc)  for  the  two  lowest  bands  in  silicon 
in  the  (111)  direction.  Even  using  the  Fermi  Golden 
Rule  rates,  there  is  a  non-zero  ionization  rate  for  states 
with  energy  less  than  1.5  eV.  For  the  collision  broadened 
case,  the  effective  threshold  is  less  than  1.2  eV.  Figure  3(a) 
shows  contours  of  equal  ionization  rate  corresponding  to 
the  average  ionization  rate  for  the  equal  energy  contours 
in  figure  3(b).  Here  kz  =  0,  and  the  states  are  in  the  first 
conduction  band.  The  wave-vector  anisotropy  is  small 
for  most  states  in  the  kz  =  0  plane.  This  is  verified  in 
figure  3;  here  the  average  ionization  rate  contours  follow 
the  energy  contours  for  almost  all  energies.  We  have 


found  deviations  of  less  than  25%  from  the  Ravc(£)  even 
for  states  with  £  <  2.0 eV.  However,  the  anisotropy 
becomes  large  below  1.5  eV,  and  may  be  play  an  import¬ 
ant  role  for  low  fields. 
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